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Abstract 

The cluster variation ~ Pade approximant method is a recently proposed tool, 
based on the extrapolation of low/high temperature results obtained with 
the cluster variation method, for the determination of critical parameters in 
Ising-like models. Here the method is applied to the three-dimensional simple 
cubic Ising model, and new results, obtained with an 18-site basic cluster, are 
reported. Other techniques for extracting non-classical critical exponents are 
also applied and their results compared with those by the cluster variation - 
Fade approximant method. 
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The cluster variation method (CVM) Hj-^j is a hierarchy of approximations which gen- 
erahzes the well-known mean field approximation and has been widely applied in the last 
decades, mainly to study the equilibrium properties of classical, discrete lattice models with 
short range interactions. The CVM results are more and more accurate as the size of the 
clusters considered increases (at least in a specific way, see below), but the critical exponents 
take always the mean field (classical) values. The issue of extracting non-classical critical ex- 
ponents from mean field approximations has been the subject of a certain amount of research 
work in recent years. As far as the CVM is concerned, a few schemes have been proposed 
in recent years to give estimates of critical exponents from the CVM results. One of these 
schemes, the cluster variation - Pade approximant method (CVPAM) |^J^ , was specifically 
devised for the CVM and exploits its great accuracy at high and low temperatures by means 
of an extrapolation of the thermodynamic quantities based on Fade approximants. 

In the present paper I report on an investigation on the CVM approximation for the Ising 
model on the simple cubic lattice with the largest basic cluster (18 sites) ever considered. 
The results of this approximation are used to give non-classical estimates of the three- 
dimensional Ising critical exponents, using mainly the cluster variation - Pade approximant 
method (CVPAM). Other schemes, like the coherent anomaly method (CAM) ||^ and an 
approach by Tome and de Oliveira [0 are also considered. 

I shall study the Ising model on the simple cubic lattice, with nearest-neighbor (NN) 
interactions only, described by the reduced Hamiltonian 

where K is the (reduced) interaction energy, Sj = ±1 is the usual Ising variable at site i, the 
summation runs on NN pairs, and ks and T are, as customary, Boltzmann's constant and 
absolute temperature, respectively. 

The CVM is a variational method based on the minimization of a free energy density 
which is obtained p|,^ by truncating the cumulant expansion of the exact variational prin- 
ciple of equilibrium statistical mechanics. The approximate free energy density depends on 
the density matrix (matrices) of the largest cluster(s) entering the expansion, which com- 
pletely determine the approximation. The simple cubic Ising model has been investigated, 
using the CVM, by Kikuchi [|l| in the NN pair, square and cube approximations and by the 
present author |0] in the star-cube approximation. Here I shall use the 18-point approxima- 
tion which is obtained by choosing as basic cluster the 3x3x2 cluster obtained by joining 
four cubes as in Fig. |l|, which is the largest cluster ever considered for this lattice. The 
choice of this cluster is motivated by Schlijper's observation [^ that, due to the existence of 
a transfer matrix for the model, accuracy can be increased by enlarging the basic clusters in 
d — \ dimensions only, where d is the lattice dimensionality. One could imagine a series (the 



generalization of the so-called C-series by Kikuchi and Brush |]I0[) of basic clusters made of 
L X L x2 sites, where for L = 2 one has the cube approximation, for L = 3 the present one, 
and for L > 3 approximations which cannot be dealt with using current computers. 

Following An [Q, the (reduced) free energy density to be minimized can be written in 
the form 

/(pis) = -3fCTr(siS2P4)+Tr(pi8lnpi8) - 2Tr(pi2 lnpi2) - Tr(p9lnp9) 

+Tr(p8lnp8) +2Tr(p6lnp6) -Tr(p4lnp4), (2) 



where pn denotes the density matrix of the n-point cluster (see Fig. |l|), while Tr(siS2P4) is 
the NN correlation (si and S2 being any NN sites on the square plaquette). 

The free energy density can be regarded as dependent on pig only, because the subcluster 
density matrices can be defined as suitable partial traces of pig. Since the hamiltonian is 
classical, the density matrices are diagonal, and hence our free energy density depends, in 
principle, on the 2^^ diagonal elements of pis corresponding to the spin configurations of 
the 18-point cluster. These elements are, however, not all independent. First of all, the 
density matrices must be normalized according to Trpig = 1. In addition, many elements 
are degenerate because of the lattice symmetries. The nondegenerate pig elements turn out 
to be 17676, a rather small number if compared with 2^^ = 262144. Taking into account 
lattice symmetries the size of the problem is reduced by a factor slightly smaller than 16, 
which is the number of elements of the symmetry group of our 18-point basic cluster. Finally, 
it must be observed that pi2 can be defined in different ways as a partial trace of pig. In 
order to ensure that these different traces yield the same density matrix one has to impose 
1134 constraints on the elements of pig. Once these constraints are satisfied no ambiguity is 
left in the definition of the other subcluster density matrices. Therefore, one is left with the 
problem of finding the minimum of a function of 17676 variables, with 1134 (leaving apart 
the trivial normalization constraint) linear constraints among them. The problem can be 
easily treated in the framework of the natural iteration method (NIM) [pl] , p!2[ , developed 
by Kikuchi for the solution of the CVM variational problem. The solution for a single value 
of K (not too close to the critical point) can be found on a modern personal computer in a 
time of the order of ten minutes. 

The accuracy of the present approximation can be assessed in several ways. For instance, 
I have compared in Tab. | the present K^ value with those from other CVM approximations 
and with the best estimates |r5|-P^. Another interesting check is the comparison, in Tab. 
|l|, of our magnetization values m = (sj) with those given by the formula by Talapov and 
Blote [O, determined on the basis of high precision simulations and finite size scaling. 



It is interesting to observe that the best agreement between the two methods occurs in 
the middle of the temperature range considered (which lies within the temperature range 
t = 1—Kc/K G (0.0005; 0.26), where Talapov and Blote regard their result as very accurate). 
For T = 1/K = 3.7 the two magnetizations differ only by 10^^. For larger temperatures our 
results are certainly less accurate than those by Talapov and Blote, while the inverse must 
be true for smaller temperatures (the result by Talapov and Blote is significantly smaller 
than 1 at very low temperatures and has a maximum around T = 1.85). 

In order to obtain non-classical estimates of the critical exponents, I shall now apply the 
CVPAM according to the rules outlined in [^Q . In the CVPAM one computes a thermody- 
namical function F for a set of temperature values in a range where the CVM approximation 
can be regarded as very accurate. These values are then used as a basis for extrapolation 
by means of Fade approximants and their generalizations. 

As a first step I shall consider the low temperature magnetization as a function of the 
variable x = e~^ . In order to determine a temperature range x < x^ax in which the 18- 
point CVM is very accurate, I compare it with a lower order CVM approximation. Although 
the star-cube approximation is slightly more accurate, I choose the cube approximation for 
this purpose, since, as I mentioned above, the cube and 18-point approximations can be 
thought of as belonging to the same series. Requiring that the magnetization difference 



is less than the empirically determined threshold e = 10~^ I obtain Xmax = 0.75. The 
function m{x) has now to be extrapolated to estimate its singular properties, taking into 
account also confluent singularities, that is corrections to scaling. To this end I use Adler's 
generalizations [[18H20[|, usually denoted by Ml and M2, of the ordinary Pade approximant 



method (I recall that an [L, M] Pade approximant is simply the ratio of two polynomials of 



degree L and M ||2l],^). Given a function F{x) with a singularity which can be assumed 
of the form (xc — x)~'*'[l + a{xc — x)^'^] method Ml considers Pade approximants to the 
logarithmic derivative of the function 

dF 

B{x) = XF{x) - {xc - x)— (3) 

for assigned Xc and A. The dominant singularity in (d/dx) lnB{x) is a pole at Xc with residue 
A — Ai if Ai < 1 and A — 1 otherwise. Method M2 considers instead, for assigned Xc and 
Ai, Pade approximants to 

G(y)^-My-l)^-, , = l-(l-_) , (4) 

which should converge to A for y = 1. In the CVPAM, the function dlnB/dx or G{y) is 
evaluated at L + M + 1 equally spaced points x„ = Xmax — n6x, x = 0,1, . . . L + M and then 
an [L, M] Pade approximant is determined by interpolation. The value of the spacing 6x 
must be adjusted empirically so that the sets of linear equations which must be solved for 
the interpolation are not badly conditioned. For the magnetization, the best conditioned 
sets of equations are obtained for 6x = 0.015. 

Applying method Ml to our magnetization estimates we have obtained the correction to 
scaling exponent Ai as a function of the critical exponent /3 for several values of the trial 



critical temperature Tc = l/Kc- It is known pO| that the plots given by Ml have a different 



curvature above and below the critical temperature. In the ordered phase the Ai versus (3 
plot is bent upward, while in the disordered phase it is bent downward. Using this criterion 
we can locate the critical temperature in the range 4.512 < T^ < 4.515, which corresponds 
to 0.22148 < Kc < 0.22163, in agreement with the most recent estimates p!3|-p!7|. From 
the corresponding plots, reported in Fig. |], one reads the estimates 0.323 < (3 < 0.332 
for the magnetization critical exponent and 0.75 < Ai < 0.82 for the correction to scaling 
exponent. As in the case of the face-centered cubic lattice ||^, the critical exponent result is 
consistent with recent estimates |T^JT^, while the correction to scaling one is substantially 
higher. Method M2 does not provide a clear-cut way to estimate T^, but using the result by 
Ml as an input we get (see Fig. |) 0.322 < (3 < 0.331 and 0.79 < Ai < 0.87. 

I tried to analyze the high temperature susceptibility data (which typically give much 
better results than the low temperature one) in the same way. I considered the susceptibility 
as a function of the variable w = tanhi^ and determined a region w < Wmax of the disor- 
dered phase in which the 18-point CVM approximation can be regarded as very accurate by 
comparing the NN correlation (which is bounded, and hence is better than the susceptibility 
for such a test) with that given by the cube approximation. The two estimates differs by less 
than e = 10~^ when w < Wmax = 0.13. Unfortunately, both methods Ml and M2 failed to 
converge and I had to resort to ordinary Pade approximants for the logarithmic derivative 
of the susceptibility. The points for the interpolation were chosen as described previously. 



with a spacing 5w = 0.003. Results for the critical exponent 7 from [L, L] approximants 



biased with the above estimates for T^ are reported in Tab. |T|. Including also results from 
[L, L ± 1] approximants I can conclude 1.237 < 7 < 1.248, which again is consistent with 
the most recent estimates. 

Finally, I have tried to analyze the magnetization data by means of other techniques, 
namely the CAM and a similar approach by Tome and de Oliveira [^. 

The CAM scaling hypothesis for the magnetization |^ is that its critical amplitude, 
defined by m{T) ~ B*{T* - Tf'"^ must diverge as B* ~ iT* - TJ^-^/^ where T^ is the 
true critical temperature and T* and B* are the estimates for the critical temperature and 
the critical amplitude in a given approximation. Using the above mentioned best estimate 
Kc = 0.22165 for the critical temperature, a fit on the results from the pair, cube, star-cube 
and 18-point CVM approximation (the plaquette approximation was discarded since it was 
clearly out of the curve) I got (3 = 0.351. Poorer results were obtained discarding, in addition 
to the plaquette approximation, the pair approximation or the star-cube approximation 
{/3 = 0.415 and 0.212 respectively). 

The approach by Tome and de Oliveira [^ is based on the scaling assumption m{Tc) ~ 
(T* — Tc)^. The best results with this approach have been found by fitting the results 
from the pair, cube and 18-point approximations {jS = 0.312) and those from the cube and 
18-point approximations only (/3 = 0.344). 

In conclusion, I have developed an 18-point (the largest maximal cluster ever consid- 
ered) CVM approximation for the simple cubic lattice and applied it to the Ising model. 
The results from this approximation have been used to extract non-classical estimates for 
the critical exponents. Among the three methods considered for this purpose, namely the 
CVPAM, the CAM and the method by Tome and de Oliveira, the CVPAM is the only 
one which gives critical exponents which are (except for the correction to scaling exponent) 
consistent with the most recent estimates. The CVPAM is therefore to be preferred when 
extrapolating CVM results to the critical region. The effort needed was essentially the same 
in all cases, since most of the labor and computer time go into the development and solving 
of the CVM approximation. In particular, the computer resources needed are remarkably 
small if compared with the requirements of extensive Monte Carlo simulations or series 
expansions, and the results are only slightly poorer, a feature which makes the CVPAM 
an interesting technique whenever powerful computers are not available and/or very high 
accuracy is not needed. 
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TABLES 



TABLE I. Comparison of various CVM estimates for the critical point with the best estimates 



Method Pair || Square |l| Cube Q Star-cube g 18-point (present) Best estimates [|T|-|T 



Kr 



0.20273 0.21693 0.21829 



0.2187 



0.2199 



0.22165 



TABLE IL Comparison of our estimate of the magnetization with that by Talapov and Blote 



l/K 



m (present work) 



m (Talapov and Blote |17|) 



3.4 
3.5 
3.6 
3.7 
3.8 
3.9 
4.0 
4.1 
4.2 
4.3 
4.4 
4.5 



0.8972562 
0.8806417 
0.8616750 
0.8399256 
0.8148173 
0.7855490 
0.7509519 
0.7092094 
0.6572414 
0.5891051 
0.4905811 
0.3067063 



0.8972440 
0.8806366 
0.8616735 
0.8399255 
0.8148161 
0.7855416 
0.7509251 
0.7091249 
0.6569722 
0.5881361 
0.4859045 
0.2378014 



L 



TABLE III. Critical exponent 7 from biased [L,L] approximants 



4.512 



4.513 



4.514 



4.515 



4 
5 
6 
7 
8 
9 
10 



1.24748 


1.24520 


1.24296 


1.24076 


1.24748 


1.24519 


1.24289 


1.24058 


1.24748 


1.24502 


1.24217 


1.23895 


1.24749 


1.24505 


1.24228 


1.23922 


1.24748 


1.24482 


1.24155 


1.23769 


1.24749 


1.24483 


1.24156 


1.23769 


1.24746 


1.24487 


1.24207 


1.23905 



FIGURES 
FIG. 1. The 18-point basic cluster and its subclusters 

FIG. 2. Ai vs p plot obtained by method Ml, for T^ = 4.512 (a), 4.513 (b), 4.514(c) and 4.515 
(d), using [L, M] approximants with 3<L<5, L — 1<M<L + 1. 

FIG. 3. p vs Ai plot obtained by method M2, for T^ = 4.512 (a), 4.513 (b), 4.514(c) and 4.515 
(d), using [L, M] approximants with 3<L<5, L — 1<M<L + 1. 
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Figure 2 
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